AD-A264  111 


Army  Research  Laboratory 


Navier-Stokes  Computations  for 
Kinetic  Energy  Projectiles  in 
Steady  Coning  Motion  - 
A  Predictive  Capability  for 
Pitch  Damping 


Paul  Weinacht 
Walter  B.  Sturek 


'^/x"  '"R'  ''X.- 

ARL-TR-112  April  1993 


APPROVED  FOR  PUBUC  RELEASE;  DISTRJBimON  IS  UNUMTTED. 


93-10488 


93  5  11  155 


NOTICES 


Destroy  this  report  when  it  is  no  longer  needed.  DO  NOT  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained  from  the  National  Technical  Information 
Service,  U.S.  Department  of  Commerce.  5285  Port  Royal  Road,  Springfield,  VA  22161. 


The  findings  of  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army 
position,  unless  so  designated  by  other  authorized  documents. 


The  use  of  trade  names  or  manufacturers’  names  in  this  report  does  not  constitute 
indorsement  of  any  commercial  product. 


REPORT  DOCUMENTATION  PAGE 


farm  Approved 
0MB  No  0704-0)89 


'eooft'Ho  Oorcjen  *or  .  oii»»ctton  nforf^.Kion  >s  r-O  4we#  tqe  ‘  '<»%OOr%e  r.i.twding  :ovtruc!i<:/'>  .vsr  '•■  ^  ;  j.j!  j  —i 

.jathe»»nq  J»'0  the  OdU  nteeded,  mo  coit'Di^tinq  .md the  a»  i«to»'?»3i*on  Send  cotntrienti  r^Odfdjrq  th,^  Durden  ,r  ir^  T»^ef  ivij-e-  t  ^ 

4.0Jle<ti0n  0^  frtTormdtJon.  rhciuding  vuqqe^tron^  tor  redutif’q  th.\  oufOen  to  ^d»h«nqton  •^eakJqudrter^  S«rv»ce\,  Direcordte  'or  htcrmdi.on  Ooerdtxrjn^  jnq  *'eoci't».  '■■*  '■*>  e?tp#Vv-r' 
Oavi^  Highway.  Su>ie  U04.  ArhnqtOh.  ■r  A  22202-4  302.  dnd  to  the  O^t'ce  O*  Menaqement  jrxi  0udqet  Pjoerworti  rteductton  Projeot  [0/04-0 ’S6>,  /^dih.nqton  vC  20SC3 


1.  AGENCY  USE  ONLY  (Leave  bUnk) 


2.  REPORT  DATE 

April  1993 


4.  TITLE  AND  SUBTITLE  S.  FUN 

NAVIER-STOKES  COMPUTATIONS  FOR  KINETIC  ENEROY  PROJECTILES 
IN  STEADY  CONING  MOTION  -  A  PREDICTIVE  CAPABILITY  FOR  PITCH 
DAMPING  PR: 

6.  AUTHOR(S) 


3.  REPORT  TYPE  AND  OATES  COVERtO 

Final,  September  1989  -  January  1992 


S.  FUNDING  NUMBERS 


1L162618AH80 
626iaA-00-001  AJ 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 

ARL-TR-112 


12b.  DISTRIBUTION  CODE 


PAUL  WEINACHT  and  WALTER  B.  STUREK 


7.  PERFORMING  ORGANIZATION  NAM£(S}  AND  AOOR£SS(ES) 

U.3.  Kcwy  Raseaxch  X.aboratozy 
ATTN:  AMSSL-NT-’PB 

Aberdeen  Proving  tSrotind,  M)  21005-5063 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 


US  Army  Research  Laboratory 
ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 

Aberdeen  Proving  Ground.  Maryland  2100S-S06G 


11.  SUPPLEMENTARY  NOTES 

This  report  supersedes  BRL-IHR-937,  March  1990, 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  releaue;  distribution  Is 
unllolted. 


13.  ABSTRACT  (Maximum  200  words) 

Previous  theoretical  investigations  have  proposed  that  the  side  force  and  moment 
acting  on  a  body  of  revolution  in  steady  coning  motion  could  be  related  to  the 

E  itch-damping  force  and  moment.  In  the  current  research  effort,  thi:'  approach  has 
een  applied  to  produce  the  first  known  Navier-Stokes  predictions  of  the  pitch 
ing  for  finned  projectiles.  The  flow  field  about  finned  kinetic  energy 
projectiles  in  steady  coning  motion  has  been  successfully  computed  using  a 
parabolized  Navier-Stokes  computational  approach.  The  computations  make  use  of  a 
rotating  coordinate  frame  in  order  to  solve  the  steady  flow  equations.  From  the 
computed  flow  field,  the  side  moment  due  to  coning  motion  is  used  to  determine  the 
pitch-damping  coefficient.  The  computational  predictions  of  the  slope  of  the  side 
moment  coefficient  with  coning  rate  normalized  by  the  sine  of  the  angle  of  attack 
have  been  compared  with  pitch  damping  coefficients  determined  from  range  firings  for 
two  finned  projectile  configurations.  The  predictions  show  good  agreement  with  the 
range  data.  This  computational  approach  provides  a  significant  predictive  capability 
for  the  design  of  kinetic  energy  projectiles  whose  terminal  ballistic  performance  can 
be  degraded  by  moderate  levels  of  yaw  at  the  target. 


14.  SUBJECT  TERMS 

Projectiles  Numerical  Analysis 

Aerodynamic  Characteristics  Computations 

Aerodynamics 
Navier-Stokes  Equations 


17.  SECURITY  CLASSIFICATION  118.  SECURITY  CLASSIFICATION  119,  SECURITY  CLASSIFICATION  I  20.  LIMITATION  OF  ABSTRACT 

OF  REPORT  I  OF  THIS  PAGE  I  OF  ABSTRACT  I 


IS.  NUMBER  OF  PAGES 

37 


16.  PRICE  CODE 


UNCLASSIFIED 


NSN  7540-01-280-5500 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  ?orm  298  :^ev  2  %9) 

Ov  4NSi  ^»q  .’39- '8 

298-^02 


Intentionally  left  blank. 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES .  v 

I.  INTRODUCTION .  1 

•2.  THEORETICAL  BACKGROUND .  3 

2.1  Relation  between  Coning  and  Pitching  Motions .  3 

2.2  Moment  Expansion .  1 

2.3  Relation  between  Side  Moment  due  to  Coning  and  Pitch  Damping  Moment.  6 

3.  COMPUTATIONAL  APPROACH .  8 

4.  RESULTS .  11 

4.1  M73.5  Projectile .  12 

4.2  MS29  Projectile .  13 

5.  CONCLUSION .  1,5 

6.  REFERENCES .  27 

LIST  OF  SYMBOLS  .  29 

Distribution  List .  31 


iooesslon  For  j  | 

RTIS 

GRA&I 

DTIC 

r.«.H 

□ 

Unuan 

□ 

Juo  1 

*  h  ;  C 

n _ 

By..---. 

Dl.'itrlt' :t  Ian/ _ __i 

Avni 

0  Codes  i 

.A- /or 

A/l 

1 

n 

Intentionally  left  blank. 


IV 


LIST  OF  FIGURES 

Figure  Page 

1  Pitching  motion  of  M735  KE  projectile  -  Round  1G123  .  17 

2  Penetration  as  a  function  of  pitch  angle,  65  degn.'e  obiir|uity  target .  17 

3  Schematic  of  coning  motion .  IS 

4  Components  of  coning  motion .  IS 

5  Schematic  of  M735  projectile  .  If) 

6  Schematic  of  M735  fin  cross-section .  19 

7  Schematic  of  M829  projectile  .  20 

8  Schematic  of  M829  fin .  20 

9  Detail  of  MS29  loading  edge .  21 

10  Detail  of  MS29  trailing  edge .  21 

11  Variation  of  side  moment  coefficient  with  coning  rate?.  M735.  Mach  1,  O)  =  2°  .  22 

12  Side  moment  coefTicient  due  to  coning,  C„^,  as  a  function  of  the  sine  of  the 

angle  of  attack,  S,  M735,  Mach  4 .  22 

Cn 

13  Development  of  normalized  side  moment  slope  due  to  coning,  over  M735 

projectile,  Mach  4 .  23 

14  Mach  number  variation  of  pitch  damping  moment  coefficient  determined  from 

side  moment  due  to  coning  compared  with  range  measurements,  M735  23 

15  Mach  number  variation  of  pitch  damping  force  coefficient  as  determined  from 

side  force  due  to  coning,  M735  .  24 

16  Fluctuating  component  of  the  swerving  motion  for  planar  pitching  motion, 

M735,  Mach  4,  one  degree  initial  yaw .  24 

17  Mach  number  variation  of  pitch  damping  moment  coefficient  determined  from 

side  moment  due  to  coning  compared  with  range?  measurenKuils,  MS29  25 

18  Mach  number  variation  of  pitch  damping  force  coefficient  as  determined  from 

side  force  due  to  coning,  MS29 .  25 


V 


Intentionally  left  blank. 


VI 


1.  INTRODUCTION 


Gun-launched  kinetic  energy  (KE)  projectiles  typically  fly  with  some  degi(‘e  <jf  pitching 
motion  caused  by  launch  disturbances  such  as  the  whipping  motion  of  the  gun  tube,  inbore 
balloting  of  the  projectile,  and  the  sabot  discard  process.  The  amplitude  of  the  pitching 
motion  decreases  (or  damps)  as  the  projectile  travels  dowurange  due  to  the  aeiodynarnic 
properties  of  the  projectile  body.  Figure  1  shows  the  pitching  motion  of  a  .NITdo  Kli  projectile 
as  observed  in  the  I3RL  Transonic  Range  and  is  representative  of  the  nuU  ion  of  fimu'd  kinetic 
energy  projectiles.  This  three-dimensional  plot  shows  the  vertical  and  horizontal  components 
of  angle  of  attack,  a  and  /3,  as  a  function  of  the  distance  dowurange.  rwo-dimeiisional 
projections  of  the  angle  of  attack  components  as  a  function  of  range  are  also  shown.  For  this 
particular  shot,  launch  disturbances  produced  initial  pitch  angles  of  more  than  si.x  degrees 
that  are  subsequently  damped  during  the  projectile’s  flight.  For  finned  KE  projectiles,  the 
rate  at  which  the  pitching  motion  is  damped  is  a  function  of  the  pitch  dam|)ing  aerodynamic 
coefficient  and  the  projectile’s  transverse  moment  of  inertia. 

The  ability  to  accurately  predict  the  pitch  damping  aerodynamic  cot'fficient  of  KE  pro¬ 
jectiles  is  of  particular  importance  to  the  projectile  designer  becau.se  the  terminal  balli.stic 
performance  of  these  projectiles  is  sensitive  to  the  pitch  angle  at  the  target.  Small  pitch  an¬ 
gles  may  result  in  significant  degradation  of  the  penetrator’s  terminal  ballistic  performance. 
Figure  2  shows  measurements  of  terminal  ballistic  penetration  as  a  function  of  pitch  angle 
for  a  long  rod  penetrator  against  a  laminated  armor  target  at  65  degrees  ol)li(iuity  (Roecker 
and  Grabarek,  1986).  Even  at  small  pitch  angles,  a  significant  loss  of  penetration  is  ob¬ 
served.  If  the  penetrator  and  the  target  are  closely  matched,  this  degradation  of  penetrator 
performance  can  result  in  the  inability  of  the  penetrator  to  defeat  the  target. 

Fin  stabilized  KE  projectiles  typically  do  not  employ  active  controls  for  damping  the 
amplitude  of  the  pitching  motion,  but  rely  instead  on  the  aerodynamic  properties  of  the 
body  to  produce  the  necessary  pitch  damping.  A  predictive  capability  for  j)itch  damping 
assists  in  the  production  of  projectile  designs  that  will  have  acceptable  le\els  of  pitch  at  the 
target.  The  development  of  such  a  predictive  capability  is  the  subject  of  this  report. 

The  pitch  damping  force  and  moments  are  generally  produced  by  the  time-dependent 
motion  of  the  body  and,  for  this  reason,  are  classified  as  dynamic  aerodynamic  derivatives. 
Despite  the  fact  that  these  coefficients  are  associated  with  the  time-de|)end('nt  motion  of  ihe 
projectile,  it  may  be  possible  to  determine  the  pitch  damping  codflcients  using  steady  mo¬ 
tion.  By  applying  linear  flight  mechanics  theory  such  a,s  that  developed  by  .Murphy  (196.3),  it 
can  be  shown  that  the  aerodynamic  side  force  and  moment  coelficients  acting  on  a  projectile 
in  steady  coning  motion  can  be  related  to  the  pitch  damping  force  and  moment  coefficients. 
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Steady  c:oiiiiig  motion  is  defined  as  tiie  motion  perfoniusl  l)y  a  missile  llx  iii!'  al  a  constant 
angle  with  respect  to  the  free  stream  vadocity  vector  aiul  uiKhagoing  a  rotation  at  a  constant 
angular  velocity  about  a  line  parallel  to  tlje  freestream  velocity  vector  and  <a/mcidenl  with  the 
projectile  center  of  gravity.  This  is  shown  schematically  in  Figure  3.  Steady  coning  motion 
can  be  decomposed  into  constant  amplitude  sinusoidal  pitching  motions  in  the  horizontal  and 
vertical  planes  (Figure  4).  Steady  coning  motion  has  the  advantag<>  of  !)('ing  a  stea<ly  motion 
when  viewed  from  llie  appropriate  coordinate  frame,  while  at  the  same  time  l)eing  composed 
of  pitching  motions.  Idie  use  of  steady  coning  motion  to  determine  the  |)itch  damping 
aerodynamic  coelficic'iits  provides  an  interesting  and  cost  effective  appiaaacli  for  determining 
the  aerodynamics  that  are  normally  associated  with  unsteady  or  time-depeiuhml  motions. 

Previously,  Tobak,  SchifF,  and  Peterson  (1969)  examined  the  aerodynamics  of  Imdies  of 
revolution  in  coning  motion  and  proposed  that  the  nondinear  aerodynamic  forces  and  mo¬ 
ments  acting  on  a  body  performing  large  amplitude  non-planar  motions  could  be  represented 
by  the  aerodynamic  forces  and  moments  produced  by  four  characteristic  motions:  ( 1)  steady 
angle  of  attack;  (2)  pitching  motion;  (3)  rolling  motion;  and  (1)  coning  motion.  Typically, 
the  linear  aerodynamic  force  and  moment  formulation  considers  only  forces  and  moments 
due  to  the  first  three  motions,  and  assumes  that  a  non-planar  motion  can  be  described  by 
the  vector  sum  of  two  independent  planar  motions.  The  addition  of  coning  motion  allows 
for  coupling  between  planar  motions  in  the  non-linear  formulation.  At  small  angles  of  attack 
where  linear  variations  of  the  aerodynamic  coefficients  are  expected,  their  non-linear  theory 
also  confirms  the  linear  theory  result  that  the  side  force  and  moment  duo  to  coning  motion 
are  related  to  the  linear  pitch  damping  coefficients. 

To  provide  additional  validation  for  the  theory,  SchifF  and  Tobak  (1970)  performed  wind 
tunnel  experiments  on  a  conical  body  undergoing  separate  or  combined  spinning  and  coning 
motions.  Their  experimental  results  showed  good  agreement  with  predictions  of  the  pitch 
damping  force  and  moment  coefficients  obtained  by  using  a  linearized  theory.  They  also 
demonstrated  that,  for  their  particular  geometry  and  flow  conditions,  the  Magnus  force 
and  moment  (variation  of  side  force  and  moment  with  spin  rate  and  angle  of  attack)  were 
negligible,  thus  the  linear  pitch  damping  coefficients  could  be  directly  determined  from  the 
side  force  and  moment  due  to  coning. 

Subsequently,  SchilT  (1972)  computed  the  supersonic  inviscid  flow  about  a  conical  body 
undergoing  coning  motion.  To  compute  the  flow  around  the  body  iii  coning  motion.  SchilT 
utilized  a  rotating  coordinate  system.  Within  the  rotating  coordinate'  frame'  the  flow  was 
steady,  thus  the  ste;ady  Euler  equations  could  be  solved.  The  governing  e-eiualious  were 
modified  to  inchide  the  centrifugal  and  Coriolis  force  terms.  The  coniputc'el  re'sults  compared 
well  with  expen'inienta!  re^siilts  and  with  c;stirnate;s  of  pitch  elamping  cex'lficie'nts  ohtaiiie'd  by 
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using  a  linear  tlieory.  Later  studies  by  Agarwal  amJ  Hakidi  (I!)78).  and  Lin  (ld7<S)  also 
employed  rotating  coonlinate  frames  to  com|)ute  the  supersonic  viscous  flow  about  conical 
bodies  in  coning  motion. 

In  this  report,  pitch  damping  predictions  for  finned  projectiles  are  oblaiiuxi  using  steady 
coning  motion.  The  supersonic  viscous  flow  field  about  these  projectiles  undergoing  coning 
motion  is  determined  computationally  using  the  parabolized  Navier-.Stokes  teclmiqiu' of  Schiff 
and  Steger  (1980).  The  computations  are  performed  in  a  rotating  coordinate  frame  similar 
to  that  employed  originally  by  SchilF  (1972).  Code  modifications  re(]uiied  to  imj)lement  the 
rotating  coordinate  frame  are  discussed.  The,se  modifications  include  the  atldilion  of  the 
centrifugal  and  Coriolis  source  terms  to  the  governing  equations  and  changes  to  the  shock 
fitting  algorithm.  From  the  computed  flow  field,  the  side  force  and  moment  due  to  coning 
motion  are  used  to  determine  the  pitch  damping  coefficients.  The  relation  between  tiie  side 
force  and  moment  due  to  coning  and  the  pitch  damping  force  and  moment  coefficients  is 
discussed  and  extensions  to  the  theory  required  for  the  analysis  of  finiu'd  projectiles  are 
described.  Results  are  presented  for  two  fielded  kinetic  energy  projectile  configurations,  the 
M735  and  MS29. 


2.  THEORETICAL  BACKGROUND 

In  this  section,  the  relation  between  coning  motion  and  pitching  motion  is  first  described. 
Following  this  description,  the  relation  between  the  side  force  and  moment  due  to  coning 
motion  and  the  pitch  damping  coefficients  is  developed. 

2.1  Relation  between  Coning  and  Pitching  Motions.  .As  was  discus.sed  previ¬ 
ously,  steady  coning  motion  is  defined  as  the  motion  performed  by  a  missile  flying  at  a 
constant  angle  with  respect  to  the  free  stream  velocity  vector  (angle  of  attack)  and  undergo¬ 
ing  a  rotation  at  a  constant  angular  velocity  about  a  line  parallel  to  the  freest  ream  velocitv 
vector  and  coincident  with  the  projectile  center  of  gravity.  This  is  shown  schematically  in 
Figure  3.  The  vertical  and  horizontal  components  of  the  angle  of  attack,  a  and  i,  vary  in  a 
periodic  fashion  as  the  projectile  rotates  about  the  free-stream  velocity  vector,  as  shown  in 
Figure  4.  The  total  angle  of  attack,  at  ^  is  constant,  however. 

Both  of  these  components  of  the  angle  of  attack,  when  plottrx!  as  a  function  of  time, 
are  sinusoidal,  constant  amplitude,  pitching  motions  that  arc  out  of  phase  with  each  other 
by  one  quarter  of  a  cycle,  as  shown  in  Figure  4.  By  dccoirqrosing  coning  motion  in  this 
fashion,  it  can  be  oijserveri  that  coning  motion  contains  a  specific  linear  combination  of  two 


orthogonal  planar  |)itcliing  motions. 

The  term  steady  coning  motion  describes  the  motion  of  the  body  about  the  freest r<'am 
velocity  vector,  but  does  not  completely  describe  the  motion  of  the  body.  In  particular,  t’u; 
projectile  may  rotate  (or  spin)  about  its  longitudinal  axis.  In  this  rejxut,  a  particular  form 
of  coning  motion,  steady  hinnr  coning  motion,  is  utilized.  In  steady  lunar  cf)ning  motifui. 
the  angular  velocity  of  the  [rrojectile  results  purely  from  the  rotation  of  llu*  [)roj(>ctile  about 
the  freestream  velocity  vector.  This  produces  a  component  of  angular  velociiy  along  the 
projectile  axis,  which,  by  definition,  is  the  spin  rate  of  th<?  projectile,  t  he  relation  betw<  ('n 
spin  rate,  p,  and  coning  rate,  4),  for  the  ca.se  of  steady  lunar  coning  motion  is  shown  Irelow. 

p  =  4>  cos  at  =  ( 1 ) 

By  specifying  both  the  coning  rate  and  the  projectile  spin  rate,  tlu;  projectile  motion  is 
now  completely  defined.  For  the  particular  case  of  steady  lunar  coning,  the  motion  can  !>e 
decomposed  into  a  combination  of  two  orthogonal  planar  pitching  motions,  plus  a  spinning 
motion  at  angle  of  attack. 

Planar  pitching  motion  is  clearly  a  tinte-dcpendent  motion  that  piodnces  a  timi'-depeiiflent 
flow  field  about  the  projectile.  Steady  lunar  cotiing  motion,  on  the  other  hand,  will  be  a 
steady  motion  when  viewed  from  a  reference  frame  that  is  attached  to  the  body  (and  there¬ 
fore  rotating  and  translating  with  the  body).  In  this  reference  system,  no  rotation  of  the 
pitch  plane  or  the  body  with  respect  to  the  reference  frame  will  be  t  l)scr\  ed.  It  is  important 
to  realize  that  the  steadiness  of  the  motion  does  not  recpiire  that  the  body  have  s])ecial  forms 
of  geometric  symmetry  (i.e.  axisymmetry).  Furthermore,  the  flow  field  generated  t)y  steady 
lunar  coning  motion  is  expected  to  be  steady,  when  viewed  from  a  coordinate  frame  attached 
to  the  body.  Steady  flow  modeling  techniques  can  be  applied  to  determine  tlx*  flow  field  due 
to  steady  lunar  coning  motion  under  the  constraints  that  both  the  coning  rate  and  the  angle 
of  attack  are  small.  (Clearly,  the  flow  may  become  unsteady  at  high  coning  rales  or  high 
angles  of  attack,  in  much  the  same  way  the  flow  over  a  body  at  fixed  angle  of  attack  at 
high  inci<lencc  can  become  unsteady  due  to  vortex  shedding.)  Since  the  reference  frame  is 
a  non-inertial  system  due  to  the  rotation  of  the  coordinate  systcun.  the  governing  equations 
must  be  modified.  Further  details  on  the  implementation  of  the  rotating  IVatne  are  provided 
in  the  di.scussion  of  the  computational  approach. 

2.2  Moment  Expansion.  It  is  common  in  aeroballistic  applications  to  utilize  a 
mi,ssile- fixed  non-rolling  coordinate  .system  to  describe  both  the  kinematics  and  the  system 
of  forces  and  moments  that  act  on  the  projectile  in  flight  (Miupliv  IDfid).  Tin'  non-roiling 
coordinate  system  affords  some  simplifications,  particularly  in  describing  the  kinematic.s. 
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Ill  this  report,  the  piiinary  reason  for  initially  deseribing  the  aei'o<lynaii lic  nuaneitts  using 
the  non-rolling  cotirdinate  system  is  the  fact  that  the  descrijitiun  is  ucll-cst ablislied.  The 
non-rolling  coordinate  frame  is  an  orthogonal  right-handed  system  {.r.i).:)  ;if<>red  at  tin* 
bo<ly  center  of  gravity.  The  is  ns<‘d  to  distinguish  the  non-rolling  axis  system  from  t  lie 
computational  coordinate  system.  The  i  axis  is  aligned  along  tiu'  projectile  longitudinal 
axis  with  the  positive  direction  oriented  towards  the  projectile  nos(>.  I'he  f  axis  is  "inii  iaily" 
oriented  downward  with  the  i  —  i  plane  perpr-ndicular  to  the  groninl.  I  he  angular  mo 
tion  of  the  non-rolling  coorditiatc'  frame  is  such  that  the  angular  ridocily  of  the  (Oordinate 
frame  with  res|)eet  to  an  inertial  frame  is  zero  in  the  direction  of  the  .i-  .ixis.  All  hough  the 
time-deperulent  orientation  of  the  non-rolling  frame  may  !><•  hard  to  visualize,  the  non-roiling 
frame  is  essentially  ecpuvalent  to  the  “fixed’’  plane  coordinate  systcmi  for  small  amplitude 
motiorns.  In  the  fixed  plane  coordinate  system,  the  x  —  z  plane  remains  perpimdiciilar  to  t  he 
ground  for  all  time.  Further  details  on  these  coordinate  frames  can  be  found  in  the  report 
by  Murphy  ( lOfid). 


The  moment  expansion  for  a  finned  missile  in  the  non-rolling  cooidiuate  fr.iine  is  shown 
in  Ecpiation  2.  This  moment  expansion  is  a  variant  of  tlie  expansion  discussed  by  Murjihy 
(1963)  for  symmetric  missiles,  'rite  most  im|K)rtant  dilh'rence  Is  that  ilie  expression  here 
includes  a  side  moment  due  to  angle  of  attack,  C„„.  The  moment  fm  imilaiion  nsi's  complex 
variables  to  separate  the  moment  components.  Cm  and  Cn,  that  arc  oriented  along  the  i) 
and  z  axes,  respectively.  The  third  moment  component,  the  roll  moment,  can  be  handled 
separately  and  is  not  of  consequence  in  this  study. 


\(  -  -\Cm,  +  iCmJ- 


(2) 


In  the  moment  expansion,  the  pitching  moment  coidficiimt ,  and  pitch  flanqnng 

moment  coefficient,  Cm^  -i-Cm^,  produce  moments  that  are  proportion, d  to  the  complex  yaw. 

and  yawing  rate  respectively.  (In  the  analysis  presented  heix'.  tln'ie  is  no  need  to 
distinguish  between  pitch  and  yaw,  and  the  terms  may  be  interchangi'd.  The  usage  fidlows 
that  of  Murphy  (1963).)  The  Magnus  moment  coelficient,  arconnts  for  a  side  moment 

due  to  flow  asymmetries  from  a  combination  of  spin  and  angle  of  attack. 

The  side  moment  due  to  angle  of  attack,  Cn.,,  is  retained  to  account  for  a  side  moment 
variation  with  angle  of  attack  that  Is  caused  by  the  beveled  fins.  'I'he  exisienre  of  this  side 
moment  was  revealed  in  a  previous  study  (Weinacht  and  Stnrt'k  1990).  'I’he  mechanism  for 
this  moment  can  be  explained  by  examining  the  wind  and  le(*-side  fins.  I’lie  fins,  wliicli  are 
beveled  to  produce  roll,  may  experience  rlifferent  flow  on  tlie  wiiul  and  lee-sides  of  the  body. 
The  lee-side  fin  may  be  immersed  in  the  wake  created  by  the  body.  If  the  How  is  snfFiciently 
different  on  the  wind  and  lee-sid(.'s,  the  bevels  on  the  wind  and  lee  fins  will  each  |)roduce 
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a  difTerent  lift  force.  The  vector  sum  of  the  forc<“s  on  these  fins  can  he  non-zeio.  llierei>\- 
producing  a  side  moment.  Of  course,  at  zero  angle  of  attack,  the  hevels  puxluce  only  a  roll 
torque;  there  is  no  side  moment  since  the  vector  sum  of  the  lift  produced  by  the  fin  hev<*!s 
is  zero. 

The  moment  expansion  presented  in  Equation  2  does  not  account  for  variations  in  tlie 
aerodynamic  coefficients  due  to  roll  orientation.  Murphy  (IDG.'J)  has  shown  that  the  form  of 
the  linear  force  and  moment  expansion  for  a  symmetric  finned  missih-  with  thia'c  (U'  more 
fins  should  have  the  same  form  as  for  a  body  of  revolution  and  that  the  elh-cts  of  roll 
orientation  are  a  higher  order  effect.  Compulations  for  finned  kinetic  energy  pnijectiles 
have  also  demonstrated  the  effect  of  roll  orientation  on  the  aerodynamic  co(dficienls  is  small 
(Weinacht  and  Sturck  1990).  In-flight  effects  of  roll  orientation  are  typically  not  observed 
because  the  projectile  spin  rate  is  greater  than  the  pitching  frerpiency,  causing  the  effects 
of  roll  orientation  to  be  averaged  out.  Flight  bodies  with  acnodynamic  coefficicnls  that 
exhibit  a  significant  dependence  on  roll  angle  may  need  to  be  treated  with  a  more  general 
aerodynamic  formulation  than  is  presented  here  (Tobak  and  Schilf  l!)7')a.  IdTob). 


2.3  Relation  between  Side  Moment  due  to  Coning  and  Pitch  Damping  Mo¬ 
ment.  In  order  to  develop  the  relation  between  the  side  moment  due  to  steady  lunar  coning 
motion  and  the  pitch  damping  moment  coefficient,  it  is  convenient  to  ix'solve  the  moment 
components  in  non-rolling  coordinates  into  moment  components  that  cause;  rotations  in  and 
out  of  the  plane  of  the  angle  of  attack.  This  relation  is  shown  below.  Here,  C’„,  is  the  in-plane 
moment  (the  moment  that  causes  rotation  of  tlie  body  in  the  plane  of  the  angle  of  attack), 
and  C„  is  the  side  moment  (the  moment  that  causes  rotations  of  the  body  out  of  the  angle 
of  attack  plane).  .Also  shown  are  relations  for  the  complex  angle  of  attack,  angular  rate,  and 
spin  rate.  Thc.se  relations,  valid  for  steady  lunar  coning  motion,  have  beem  simplified  from 
the  general  ca.se  of  arbitrary  motion  (Levy  and  Tobak  1970). 


Cm  +  iCn  — 


d(']) 


70 


(3) 


The  moment  formulation  cast  in  terms  of  the  in-planc  and  side  monumts  can  be  written  as 
follows. 


6'„.  +  i.c\  =  Cmj  +  +  c^j  + 


i-q  /'  ^  O  i  J 


f) 


As  expcctotl,  tlio  resulting  expression  for  the  in-pl<vn<‘  ami  side  moincnls  is  in<i<'p<Mident 
of  time.  The  in-[)lane  moment  results  only  from  the  pitching  moment,  wiiih?  tli<,‘  tcjtal  side 
moment  consists  of  contributions  from  the  side  moment  due  to  nngh'  of  attack,  Magnus 
moment,  and  pitch  damping  moment. 


Utilizing  Equation  4,  the  variation  of  side  moment  witli  coning  rate  can  iio  obtained. 
This  relation  is  valid  for  linear  variations  of  side  moment  with  coning  rate. 


a.  s 


dCn  -  Cr^J 


a(f) 


iL 

V 


—  6{'iCnp^  +  +  lUr,,,,]) 


i^] 


The  relation  of  the  side  moment  due  to  coning,  to  the  pitch  tlaini)ing  moment  coef¬ 
ficient,  [Cm,  +  Magnus  moment  coefficient,  Cn^^,  is  similar  to  that  presented 

by  Schiff  and  Tobak  (1970)  for  bodies  of  revolution.  However,  for  the  case  of  the  finned  pro¬ 
jectiles  with  beveled  or  canted  fins,  the  evaluation  of  C„  requires  tliat  the  net  side  moment. 
C„,  be  determined  at  two  separate  coning  rates  (which  may  include  zero  coning  rate),  due 
to  the  presence  of  a  side  moment  at  zero  coning  rale,  C„^(5.  For  IkxHc.’s  of  la.’volut  ion.  the 
side  moment  at  a  single  non-zero  coning  rate  is  sufficient  to  determiiK;  tlie  slope.  C„  .  since 
the  side  moment  at  zero  coning  rate  is  zero. 


Unfortunately,  the  pitch  damping  moment  coefficient  is  not  directly  related  to  the  side 
moment  due  to  coning,  because  of  the  presence  of  the  Magnus  moment  term.  Further¬ 
more,  direct  evaluation  of  the  Magnus  moment  us'’^g  a  Navicr-Stokes  approach  would  be 
difficult  and  expensive  because  an  unsteady  time-accurate  calculation  is  required  for  non- 
axisymmetric  bodies.  However,  the  Magnus  moment  coefficient  is  ty])ically  much  smaller 
than  the  pitch  damj^ing  coefficient  for  many  projectiles.  This  has  been  confirmefl  by  apply¬ 
ing  simple  inviscid  theories  to  estimate  the  Magnus  moment  and  pitch  damping  coefficients 
for  finned  projectiles  similar  to  those  examined  in  this  study  (Devan  I9S9).  Additional  con¬ 
firmation  can  be  found  from  ballistic  range  testing  of  kinetic  energy  projectiles.  The  Magnus 
moment  coefficient  can  be  quite  difficult  to  measure  due  in  part  to  its  small  magnitude 
in  relation  to  the  other  aerodynamic  coefficients  and  due  to  the  low  s|)in  rates  that  these 
projectiles  experience  in  flight. 


For  the  case  where  the  pitch  damping  moment  coefficient  is  much  larger  than  the  Magnus 
coefficient.  Equation  5  can  be  simplified  without  a  significant  loss  of  accuracy,  and  the 
following  relation  is  obtained. 

+  (fi) 


Because  this  expre.ssion  is  only  valid  in  the  linear  aerodynamics  regime  (small  angles  of 
attack),  the  cosine  of  the  angle  of  attack,  7,  is  approximately  erjual  to  otu'  and  no  longer 
appears  in  the  equation. 
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A  similar  ex[)ressi<)n  relating  side  force  due  to  coning  to  tlie  pitch  damping  foice  and 
Magnus  force  can  be  developed  using  the  same  ap|)roacli  as  discussed  al)ove.  Since  the 
Magnus  force  coelRcient  is  also  much  smaller  than  the  pitch  damping  force  coc'llicieni.  it  can 
be  ignored,  and  the  pitch  damping  force  coefficient,  C',v,  +  C'.v^  can  be  directly  odated  to  the 
side  force  coefficient  due  to  coning,  Cy^fS. 

[Cn,  +  C',v„]  ~  (7) 

Equations  6  and  7  form  the  basis  for  determining  the  pitch  damping  force  and  moment 
coefficients  in  this  research  effort.  The  side  force  and  moment  acting  on  the  [)roJect  ih?  due 
to  steady  lunar  coning  motion  can  be  determined  by  computing  the  flow  field  around  the 
projectile  and  then  integrating  the  pressure  and  shear  stres.ses  acting  on  the  body.  Using 
these  two  equations,  the  pitch  damping  coefficients  are  obtained  from  the  side  force  and 
moment. 


3.  COMPUTATIONAL  APPROACH 


Computation  of  the  viscous  flow  field  about  the  finned  projectile  configurations  was 
accomplished  by  solving  the  thin-layer  Navier-Stokes  equations  using  a  ivarabolized  .N’avier- 
Stokes  technique.  The  computations  have  been  performed  in  a  rotating  coordinate  frame 
that  is  attached  to  the  projectile  body  and  rotates  at  the  coning  rate  of  tlie  proje'ctile.  The 
fluid  flow  relative  to  the  rotating  coordinate  frame  does  not  vary  with  time,  allowing  the 
steady  (non-time  varying)  Navier-Stokes  equations  to  be  applied.  The  solution  of  the  steady 
Navier-Stokes  equations  can  be  performed  at  a  reasonable  computational  cost,  typically  one 
and  a  half  to  two  hours  on  a  Cray-2  or  Cray  X-MP  computer.  In  order  to  implement  the 
rotating  coordinate  frame,  the  governing  equations  have  been  modified  to  include  the  effect 
of  centrifugal  and  Coriolis  forces.  The  steady  thin-layer  Navier-Stokes  equations  are  shown 
below. 


M  ^  ^  //  _ 

(94  di)  ^  (X  ~  fie  (X 


iS) 


Here,  E,  F,  and  G'  are  the  inviscid  flux  vectors,  S  is  the  viscous  flux  vector,  and  II  is  the 
source  term  containing  the  Coriolis  and  centrifugal  force  terms  that  ixjsult  from  the  rotating 
coordinate  frame.  Each  of  these  matrices  are  functions  of  the  dependent  variables  reprc’sented 
by  the  vector  pn,  pv,  pw,e),  where  p  and  e  are  the  density  and  the  total  energy  per  unit 
volume,  and  u,  v,  and  tv,  are  the  velocity  components  in  tiie  x,  y,  and  r  directions.  The  flux 
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terms  are  defined  as  follows. 


pU 

pV 

pu- 

puU  +  Gp 

puV  +  Tjj.p 

pnir  +  Gp 

^=7 

pvU 

^  =  7 

PVV  +  f/gP 

G=J 

pell-'  +  Gl> 

pwU 

pioV  +  r/,p 

pndl'  +  G'P 

.  i^  +  pW  , 

r 

(e  +  p)V' 

(c  +  p)!!' 

II  =  \ 


0 

Ih 

Ih 

lU 

ih 


S’  =  -!- 


0 
dll 

+  '»2G 
(Id 

m,— +  7»2(„ 

dw 

^«l-^  +  mzG 

m3 


(fij 


where 


//j  =  — 2Qc  si'i  o-t/JD  — sin^  a<(x  —  Xcg)  + sin  Oj  COSO, 

1/3  =  ‘lOcSinoj/Di  -  2ncCOsa(/>iy  - /)n^ysin^O(  -  pQ^Z/cos^n, 

7/4  =  20c  cos  a(/)y  + /’fig  sin  Q;(  cosa,(i  -  leg)  - cos*  Q, 

7/5  =  (  —  0^  sin^  a((x  —  Xgg)  +  O^x  sin  a<  cos  ai)pu  —  (0^//siri^  Oj  +  0^!/  cos*  a,  )pv 

+(0^  sinofj  cos  at{x  -  icg)  -  O^zcos^  oit)pw  (10) 


U  =  t(G 

V  =  u77j.  4-  vpy  +  wp. 

W  =  uG  +  dG  +  "-’G 

(H) 

rrii  =  (p  +  P()(Cr +Cy  +G^) 

1  .  X,  >  ^  dv  .  city  ^ 

m,  =  j(,,  +  ,„)(C,_  +  C, 

^  /  P  Pt  v/,3  .2  1  di\‘ 

+m2(uSx  +  eCy  +  luG) 

(1-2) 

^2  _  TP 
p 

(13) 

4-  D^  + 

(M) 

G  = 

r}z  =  -  y(.H)  Vy  =  P.  =  -/(  -Xi/yc) 

G  “  ~  X^/lrj)  Cy  ~  •7(~X^Cg)  G  “  '7(X(;Z/g) 

•7  =  l/(x^{?/gZ(;  -  p<;Zg)) 

(15) 
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The  form  of  tlie  source  terms,  ll y.  Hi,  If.x,  II5,  assumes  that  the  x-axis  is  orieutc'd  along 
the  projectile's  longitudinal  axis  and  tlie  x  —  z  plane  is  in  the  pitch-plane. 

The  pressure,  p,  can  be  relateil  to  the  dependent  variables  by  ajiplying  tin?  ideal  gas  law. 

p  =  (7  -  1)  e-  ( IG) 

The  turbulent  viscosity,  /it,  which  appears  in  the  viscous  matrices,  is  ccun[)uted  using  the 
turbulence  model  of  Baldwin  and  Lomax  (I97S). 

The  thin-layer  ccpiations  are  solved  using  the  paraboliztxl  Xavu-r-Stokes  l<-clnii(iue  of 
SchifT  and  Steger  (lOSO).  Following  the  approach  of  Schilf  aiul  Sieger,  the  governing  ecpia- 
tions,  which  have  been  modified  here  to  include  the  Coriolis  and  centrifugal  force  terms, 
are  solved  using  a  conservative,  approximately  factored,  implicit  finite-dilTercnce  numerical 
algorithm  as  formulated  by  Beam  and  Warming  (1978).  The  equations  are  first  linearized 
and  placed  in  delta  form,  where  tlie  equations  are  solved  for  the  diirerence  in  the  dependent 
variables  rather  than  the  variables  themselves.  This  set  of  eijuations  is  tlu'ii  factored  using 
the  approach  of  Beam  and  Warming.  The  following  set  of  equations  is  obtained. 

/i;  +  (i  +  =  fins  (i7) 

+  ( l  -  o)A?  (ifCi  -  ^  {?(/»'))]  ( I  s) 

RHS  =  -(/ii  -  Ai-'W  +  a{E’,  -  -  [({././)'+' -  iiJ.n’Ei'' 

+<<  [(g‘(£/-/)'  +  Q*'(F/Jy  +  cr'(c/./v] 

+  H' - +<t)  (I!)) 

lie 

The  form  of  the  equations,  as  well  as  the  notation,  is  similar  to  that  used  by  SchifT  and 
Steger.  Here,  A,  D,  C,  and  M  are  the  .Jacobian  matrices  of  the  flux  vectors  F.  /',  G. 
and  S.  Further  details  on  the  definitions  of  these  matrices  can  be  found  in  the  pajier  of 
SchifF  and  Steger  (1980).  The  important  dilTercnce  in  the  current  formulation  is  th<’  aridition 
of  the  matrices  D  and  ll  due  to  the  rotating  coordinate  system.  .411  hough  the  .lacobian 
matrix  of  the  source  term,  D,  can  be  included  in  either  the  circumferential  inversion  or  in 
the  normal  inversion,  including  this  term  in  the  circumferential  inversion  simplifies  slightly 
the  implementation  of  the  shock  fitting  boundary  conditions. 

The  computations  presented  here  were  performed  using  a  shock  fitting  procedure  re¬ 
ported  by  Rai  and  Chaussee  (198.3).  This  procedure  solvx's  the  fi\'e  Rankiiu'-llugoniot  jump 
conditions,  two  geometric  shock-propagation  conditions,  and  ou(’  compatibility  eipiation  to 
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determine  tlie  valiu's  of  the  five  de[)<Mu!ent  variables  iimncdiately  downsi  ream  (jf  the  sliork. 
as  well  as  the  position  of  the  sliock.  By  including  the  implicit  part  of  tlie  sonrte  term  due 
to  the  rotating  coordinate  frame  in  the  circumferential  inversion,  the  shock  lilting  procedure 
of  Rai  and  Chaussec  can  he  used  without  modification  as  long  as  the  c(arecl  free  stream 
conditions  are  specified.  The  freestream  values  of  the  (h'pendent  variables  are  show  n  irelow 
in  non-dimensional  form. 

P  ~  ' 

pu  =  cos  a  I  ijQ^.  sin  Oj 

pv  =  llel-  at  —  (x  —  Xcg)  sin  «() 
pw  =  i\/oo  sin  Oi  —  yfic  cosa^ 

e  =  Px)/(7  -  0  +  ^{(^'^ooCOSOt  +  }/f)cSina,)^ 

+(Q<-(- cos  Oi  —  (j:  —  Xcg)s\i\  at)f  +  (i'l/.^sinO(  —  cos  r>, )‘ }  (20) 

The  computational  results  presented  here  were  obtaimal  using  a  grid  that  consisted  of 
60  points  between  the  body  and  the  bow  shock.  Due  to  a  lack  of  circumferential  symmetry, 
gridding  was  performed  around  the  entire  circumference  of  the  body.  Over  llu’  forebody,  72 
circumferential  points  were  used.  Grid  resolution  was  increased  to  200  points  on  the  fin  hub. 
The  grid  over  the  finned  portion  of  the  body  was  generated  using  an  elliptic  grid  generation 
scheme  presented  by  Rai,  Chaussee,  and  Rizk  (1983).  On  the  a.xisyinm(d,ric  part  of  the  body, 
about  50  marching  (a.xial)  planes  were  required  for  each  caliber  of  body  length.  A.xial  grid 
resolution  was  doubled  over  the  finned  portion  of  the  body. 

4.  RESULTS 

Computations  have  been  performed  to  determine  the  aerodynamics  of  kinetic  energy 
projectiles  in  steady  lunar  coning  motion.  Results  have  been  obtained  for  two  fielded  kinetic 
energy  projectiles;  the  .M735  and  the  M820.  Schematics  of  these  project  ik's,  including  rletails 
of  the  fin  geometry,  are  shown  in  Figures  5-10.  The  fins  on  both  of  thes('  project  ih's  have 
roll-producing  beveled  surfaces  at  the  leading  and  trailing  edges  of  the  fins.  Particular  care 
has  been  taken  to  model  the  fin  geometry  accurately.  It  should  lx;  not('d  that  the  fins  on  the 
M829  projectile  overhang  the  base.  This  aspect  of  the  projectile  u-as  moileled  by  extending 
the  base  so  that  it  was  aligned  with  the  trailing  edge  of  the  fin  blades,  d'his  allowcvi  the  flow 
field  to  be  computed  up  to  the  trailing  edge  of  the  fins.  Howawer,  when  tlx-  pr«'ssure  and 
viscous  stres.ses  were  integrated  to  comptite  the  forces  acting  on  the  body,  the  contribution 
from  this  part  of  tlu;  body  was  not  considered.  Because  the  flow  is  supersonic  and  the  fins 
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are  not  immersed  in  the  recirculating  flow  in  the  l)ase,  the  flow  field  adjacent  to  this  r<‘gion 
can  be  considered  to  be  reasonably  well  nuKlole<l.  Thougli  not  shown,  the  cylindrical  portion 
of  these  bodies  have  a  number  of  circumferential  grooves  which  cover  nearly  two-thirds  of 
the  body.  The  effect  of  these  grooves  is  not  modeled  in  the  current  computations,  though  it 
is  a  subject  of  current  research. 

The  computations  have  been  performed  over  a  range  of  Mach  numl)ers  (.1/^  =  'bO  to  5.5), 
coning  rates  (oDfV  —  0.0  to  0.010),  and  angles  of  attack  (ot  =  1°  to  5°)  for  free-llighi,  sea- 
level  atmospheric  conditions.  The  variation  of  the  side  force  and  moment  with  coning  rate 
normalized  by  the  angle  of  attack  has  been  used  to  determine  the  pitch  dam[)ing  coefficients 
for  these  finned  projectiles.  Comparisons  are  made  with  data  obtained  from  range  firings 
(Brandon).  The  results  for  each  projectile  are  discussed  separately  in  the  following  two 
sections. 


4.1  M735  Projectile.  The  computed  variation  of  the  sirlc  monuMit  coefficient  with 

coning  rate  at  .Mach  1  and  two  degrees  angle  of  attack  is  shown  in  Figure  1  1.  The  variation 
of  the  side  moment  coefneient  with  coning  rate  is  seen  to  be  linear  across  the  range  of 
coning  rates  examined  here.  This  range  of  coning  rates  is  representative  of  the  pitching 
frequencies  experienced  by  the  M735  projectile  in  flight.  At  Mach  -1,  the  non-dimensional 
pitching  frequency  of  the  projectile  is  0.004,  where  the  form  of  the  non-dimensiona!ization  is 
the  same  as  for  the  coning  rate.  The  results  also  show  the  existence  of  a  small  non-zero  side 
moment  coefficient  at  zero  coning  rate.  .As  discussed  previously,  this  side  moment  is  due  to 
bevels  on  the  fins  (Weinacht  and  Sturek  1990).  The  existence  of  this  side  inoment  at  zero 
coning  rate  requires  that  computations  be  performed  for  at  least  two  coning  rates  in  order 
to  evaluate  the  variation  of  the  side  moment  coefficient  with  coning  rate.  C,,^. 

Figure  12  shows  as  a  function  of  S  (the  sine  of  the  angle  of  attack)  at  .Mach  4.  The 
dashed  line  displayed  on  this  figure  is  representative  of  a  linear  variation  of  with  6  across 
the  range  of  angles  of  attack  examined.  The  computed  results  show  that,  at  small  angles  of 
attack,  Cn^  varies  linearly  with  S,  but  depart.s  from  a  linear  variation  as  the  angle  of  attack 
increases. 

Figure  13  shows  the  development  of  Cn^l^  over  the  M735  kinetic  energy  [)rojectile  at 
Mach  4  and  two  degrees  angle  of  attack.  As  discussed  previously,  should  be  a  reason¬ 

able  representation  of  the  pitch  flamping  coefficient,  Cm^  +  Cm^,  in  the  linear  aerodynamic 
regime.  This  figure  shows  that  the  fins  contribute  most  of  the  side  monuuit  due  to  coning 
(and  hence,  the  pitch  damping)  with  a  smaller  contribution  from  the  nos*'. 

The  Mach  variation  of  C’m,  +  Cm^  ihe  M73.5,  as  detfu-mined  from  is  shown 
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in  Figure  14.  The  computed  results  are  compared  with  raiigt;  nu'asurenKMits  of  tlie  pitch 
damping  coelficient.  'I'liough  tlie  range  data  sliown  here  are  coiisidfued  weil-delerinined, 
some  scatter  is  still  evident  because  damping  rates  are  typically  dilhcult  to  m(*asure.  The 
experimental  results  do  reflect  the  expected  level  of  accuracy  in  deteriniiiing  this  coeflicient 
experimentally.  The  comparisons  show  that  the  computational  results  are  within  tlie  ac¬ 
curacy  of  the  experimental  data  and  provide  a  measure  of  validation  of  the  computational 
approach. 

The  predicted  v'ariation  of  the  damping  force,  C’.v,  +  ‘"is  function  of  .Mach  num¬ 
ber  is  shown  in  Figure  15.  This  coefficient  was  determined  from  the  variation  of  the  side 
force  coefficient  with  coning  rate,  The  pitch  damping  force  coefficient  appears  in  the 

swerve  equation  (the  equation  that  describes  in-flight  motion  of  the  projectile  center  of  grav¬ 
ity).  For  finned  projectiles,  the  fluctuating  component  of  the  swerving  motion  is  composed 
of  contributions  from  the  lift  and  pitch  damping  force.  The  fluctuating  component  of  the 
swerving  motion  for  a  damped  planar  pitching  motion  is  shown  in  Figure  16.  The  relative 
contributions  from  the  lift  and  pitch  damping  forces  are  also  shown.  Tlie  contriimtion  from 
the  pitch  damping  force  coefficient  is  seen  to  be  a  relatively  small  portion  of  the  motion, 
despite  the  large  magnitude  of  the  coefficient  itself.  Because  of  this,  it  is  very  poorly  deter¬ 
mined  from  range  firings,  thus  no  experimental  data  is  shown.  This  coefficient  is,  however, 
required  for  determining  the  change  in  the  pitch  damping  moment  coefficient  due  to  changes 
in  center  of  gravity  location.  Thus,  the  ability  to  predi'  the  pitch  damping  force  coefficient 
is  still  a  significant  result. 

4.2  M829  Projectile.  S  imilar  computations  were  performed  for  the  ,MS2y  kinetic 

energy  projectile.  Figure  17  shows  the  predicted  variation  of  the  pitch  damping  moment 
coefficient  with  .Mach  number  for  the  M829.  Again,  the  computed  results  are  compared  with 
range  measurements.  The  range  data  has  considerable  scatter  because  the  total  angle  of 
attack  experienced  by  the  projectiles  during  the  range  firings  was  very  small  (typically  less 
than  one  degree).  Thus,  the  rate  at  which  the  amplitude  of  the  pitching  motion  decreased  in 
flight  was  difficult  to  determine.  The  computational  results  are  within  the  scatter  of  the  range 
data.  Both  the  computations  and  the  range  results  show  the  order  of  magnitude  increase  in 
the  coefficient  comparerl  with  the  predictions  shown  for  the  M7.35.  This  increase  is  primarily 
due  to  the  larger  length-to-diametcr  ratio  (L/D)  of  the  MS2'J  (L/D  ~  23)  compared  with 
the  M7.35  (L/D  %  14).  The  predicted  variation  of  the  pitch  dami^ing  force  coefficient  as  a 
function  of  Mach  number  is  shown  in  Figure  IS.  Again,  no  range  data  are  shown  because 
this  coefficient  is  poorly  determined.  This  coelficient  is,  however,  useful  for  determining  the 
variation  of  the  i)itch  damping  moment  coefficicuit  with  changes  in  c<'nl<.T  of  gravity  position. 
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One  situation  wlion'  changes  in  the  center  of  gravity  location  are  often  (“nconnleretl  is 
in  aerodynamic  range  testing  of  kinetic  energy  projectiles.  Because  of  restrictions  on  firing 
kinetic  energy  projectiles  with  heavy  metal  penetrators  through  aerodynamic  t('st  ranges, 
surrogate  projectiles  are  often  used.  Externally,  these  projectiles  apitear  the  same  as  the 
fielded  round,  though  the  internal  heavy  metal  core  has  been  replaced  with  a  steel  core. 
Replacing  the  penetrator  core  can  result  in  a  shift  in  the  center  of  gravity.  In  the  case  of  the 
MS'29,  this  shift  is  greatt'r  than  a  quarter  of  a  caliber. 

The  effect  of  center  of  gravity  (CG)  shift  on  the  <lamping  coefficients  was  investigated 
in  two  ways.  First,  flow  field  computations  were  performed  to  determiiu'd  the  elfect  of  GG 
location  on  the  side  force  and  moment  due  to  coning  (and  hence  on  the  pitch  damping  force 
and  moment  coefficients).  Since  the  projectile  rotates  about  the  CG  position,  the  approach 
involved  separate  flow  field  calculations  for  each  CG  position.  The  CG  position  was  moved 
fore  and  aft  of  the  ba^seline  CG  position  by  1  and  2  body  diameters.  These  results  are 
displayed  in  Table  1  for  .Mach  numbers  of  4,  and  5.  The  pitch  damping  moment  becomes 
more  negative  (indicative  of  stronger  damping)  as  the  CG  position  move's  forward. 

Table  1.  CG  Translation  -  Comparison  of  Results  Using  CG  Translation  Faiuations  and 
Direct  Computation 


MACH 

NUMBER 

^cg 

COMPUTED 

TRANSLATED 

msm 

■aRH 

3.0 

-2.0 

224.37 

-1550.1 

224.61 

-1551.2 

-1.0 

241.89 

-1853.2 

242.01 

-1853.9 

276.93 

-2564.6 

276.81 

-2.563.7 

294.45 

-2972.8 

29'1.21 

-2970.8 

4.0 

-2.0 

179.59 

-1373.4 

179.85 

-1374.9 

-1.0 

195.32 

-1615.7 

195.45 

-1616.6 

+  1.0 

226.78 

-2194.7 

226.65 

-2193.6 

+2.0 

242.52 

-2531.5 

242.25 

5.0 

-2.0 

136.80 

-1177.2 

137.08 

-1.0 

150.78 

-1361.1 

150.92 

■m 

-1812.7 

178.62 

+2.0 

192.75 

-2080.6 

192.46 

As  an  alternative  to  the  first  approach,  the  effect  of  CG  position  on  the  j)itch  damping 
force  and  moment  coefficients  can  be  determined  using  the  center  of  gravity  translation 
relations  preisented  by  .Murphv  (1963).  The  relations,  present<'d  by  Murphy  for  tlu'  individual 


aerodynamic  coidncient.s,  are  coml)iiied  to  ol)tain  ndations  for  t  he  (xx'fficieiit  s  of  interest  liere>. 
These  relations  for  the  pitch  damping  and  Magntis  force  and  monuMits  are  siiou  n  Ijelow.  The 
coefficients  on  the  left-liand  side  of  the  equations  (denoted  with  a  “*")  re|)r<'sent  lln*  predicted 
value  for  a  CG  shift  of  Scg  body  diameters.  The  aerodynamic  cordficienls  on  tiu'  right  side 
of  these  expressions  reiuesent  the  values  for  the  bcaseline  configuration. 


—  C|V,  +  C,v.,  +  -s 

+  Gm,i 

=  C,n^  -I-  C'„i„  —  .' 

>c;i{Csg  +  CV„)  + 

= 

Since  the  pitch  damping  force  and  moment  coefficients  are  related  to  the  side  force  and 
moment  due  to  coning,  (Equations  6  and  7),  the  side  force  and  moment  coefficients  due  to 
coning  exhibit  the  same  variation  with  CG  position. 


Cy^/6  =  Cs'^/6  + 


(22) 


The  derivation  of  ihe.se  relations  does  not  require  that  the  Magnus  force  and  moment  be 
neglected.  The  effect  of  the  Magnus  coefficients  is  simply  absorbed  into  the  coelficients  Cy^ 
and  Cn^,  and  adds  no  additional  terms  to  the  right-hand  side  of  these  eciuations. 

Using  these  rotations  and  the  aerodynamic  coefficient  predictions  for  the  baseline  config¬ 
uration,  predictions  of  the  side  force  and  moment  variation  due  to  coning  for  varying  CG 
position  were  obtained.  These  results  are  included  in  Table  1.  The  dilfinences  between  the 
direct  computation  of  the  side  moment  coefficient  at  the  various  CG  locations  and  the  values 
obtained  from  the  CG  translation  relations  is  less  tlian  0.2%,  and  jirovides  additional  valida¬ 
tion  of  the  computational  approach.  The  side  moment  (and  hence  the  pitch  damping)  varies 
by  more  than  ±  3.')%  for  a  center  of  gravity  shift  of  ±2  calibers.  Changing  the  CG  position 
is  one  possible  approach  for  increasing  or  decreasing  the  pitch  damping  of  these  |uojectiles. 


5.  CONCLUSION 

The  flow  field  about  finned  kinetic  energy  projectiles  in  steady  coning  motion  has  been 
successfully  comptited  using  a  parabolized  Navier-Stokes  computational  approach.  The  com¬ 
putations  make  use  of  a  rotating  coordinate  frame.  Relative  to  this  coordinate;  frame,  the 
flow  does  not  vary  with  time,  allowing  the  steady  flow  equations  to  be  solved.  Using  lin¬ 
ear  flight  meclianics  theory,  the  side  moment  due  to  coning  is  related  to  the'  pitch  damping 
and  Magnus  mormnit  coefficients.  For  small  Magnus  moment  coc'fficie'iits.  the  pitch  damping 
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coefficit'nl  is  diroctly  proportional  to  tlie  side  moment  <lne  to  coninj!,.  I’lie  <  oini)Utational 
results  show  that  the  side  moment  due  to  coning  varies  linearly  svith  coning  rate  over  a 
range  of  coning  rates  that  encompasses  the  pitching  frequency  of  the  projectile.  For  con¬ 
stant  coning  rate,  the  side  moment  coefficient  al-so  shows  a  linear  behavior  with  the  sine  of 
the  angle  of  attack  up  to  about  two  degrees  angle  of  attack,  and  deviates  slightly  from  the 
linear  behavior  at  higher  angles  of  attack.  The  computational  pr<'dictions  of  tin-  slo[j(?  of 
the  side  moment  coefficient  with  coning  rate  normalized  by  the  sine  of  the  aiigle  of  attack 
have  been  compared  with  pitch  damping  coefficients  determined  from  range  firings.  For  the 
M735  kinetic  energy  projectile,  the  predictions  are  in  good  agreement  with  the  range  data. 
The  computational  results  for  the  MS‘29  are  within  the  scatter  of  the  range  data.  Both  the 
computational  predictions  and  the  range  data  for  the  MS29  show  a  substantial  increase  in 
the  pitch  damping  coefficient  when  compared  to  the  damping  of  the  M735.  This  is  primarily 
due  to  the  larger  length-to-diameter  ratio  of  the  M829. 

The  favorable  compari.sons  with  range  data  and  the  efficiency  of  the  computational  ap¬ 
proach  demonstrate  the  utility  of  this  newly  developed  capability.  To  date,  several  advanced 
kinetic  energy  projectile  concepts  have  been  examined  using  this  capability. 
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Figure  5.  Schematic  of  M735  projectile. 
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Figure  6.  Schema  tie  of  M735  fin  cross-,5eclioa. 
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ALL  DIMENSIONS  IN  CALIBERS  (ONE  CALIBER  =  27,05  mm) 


ALL  DIMENSIONS  IN  CALIBERS  (ONE  CALIBER  =•  27,05  mm) 
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Figure  10.  Detail  of  MS29  trailing  edge. 
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Figure  11.  Variation  of  side  moment  coefficient  with  coning  rate,  M735,  Mach  4,  q,  =  2°. 
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Figure  12.  Side  moment  coefficient  due  to  coning,  as  a  function  of  the  sine  of  the 
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Figure  13.  Development  of  normalized  side  moment  slope  due  to  coning,  —f-,  over  M735 
projectile,  Mach  4. 
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Figure  14.  Mach  number  variation  of  pitch  damping  moment  coefTicient  determined  from 
side  moment  due  to  coning  compared  with  range  measurements,  M735. 
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side  force  due  to  coning,  M735. 


Figure  16.  Fluctuating  component  of  the  swerving  motion  for  planar  pitching  motion, 
M735,  Mach  4,  one  degree  initial  yaw. 
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Figure  17.  Mach  number  variation  of  pitch  damping  moment  coeiricicnt  dcUniniiicd  from 
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Figure  18.  Mach  ruimber  variation  of  pitch  damping  force  coefficient  as  (lct<'riiiiii('(!  Trorn 
side  force  due  to  coning,  M829. 
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speed  of  sound 

pitching  moment  coefficient 

slope  of  the  pitching  moment  coefficient  with  angle  of  attack 
pitch  damping  moment  coefficient 
side  moment  coefficient 

slope  of  the  side  moment  coefficient  with  angle  of  attack 
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Magnus  moment  coefficient 

slope  of  the  normal  force  coefficient  with  angle  of  attack 

pitch  damping  force  coefficient 

slope  of  the  side  force  coefficient  with  coning  rate 

Magnus  force  coefficient 

projectile  diameter 

Jacobian  matrix  of  the  source  term,  // 

total  energy  per  unit  volume,  non-dimensionalized  by 

flux  vectors  in  transformed  coordinates 

source  term  resulting  from  rotating  coordinate  frame 

Jacobian 

characteristic  length,  typically  the  projectile  diameter 
freestream  Mach  number 

pressure,  as  used  in  thin-layer  Navier-Stokes  equations, 
non-dimensionalized  by  Poof^l^ 

spin  rate,  as  used  in  roll  equations  and  aerodynamic  coefficients 
Prandtl  number 
turbulent  Prandtl  number 

vector  of  dependent  variables  in  Navier-Stokes  equations 

total  velocity  of  fluid 

Reynolds  number,  a^oPooDJ 

distance  downrange 

center  of  gravity  shift,  calibers 

viscous  flux  vector  in  transformed  coordinates 

reference  cross  sectional  area  of  projectile, 

time 

velocity  components  in  x,  y,  and  c  directions, 
non-dimensionalized  by  a^o 

Coiitravariant  velocities  of  the  transformed  Navier-Stokes  e(|uations 
freestream  velocity  used  to  non-dimensionalize  the  s(hu  rate  and 
the  aerodynamic  coefficients 

Cartesian  coordinates  with  respect  to  the  body,  non-dimensionalized  by  D 
axial  location  of  projectile  center  of  gravity  with  respect  to  Cartesian 
coordinate  system 
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Greek  Symbols 


Of 

at 

7 

7 

<5 

Ih 

P 

4>D 

V 


vertical  component  of  total  angle  of  attack  in  non-rolling  coordinates 
total  angle  of  attack, 

horizontal  component  of  total  angle  of  attack  m  non-rolling  coordinates 

ratio  of  specific  heats,  as  used  in  Navier-Stokes  equations 

cosine  of  the  angle  of  attack,  as  used  in  aerodynamic  force 

and  moment  formulations 

sine  of  the  total  angle  of  attack 

laminar  viscosity 

elFective  turbulent  viscosity 

transformed  coordinates  in  Navier-Stokes  equations 

complex  quantity  representing  the  components  of  the  sine  of  the  angle 

of  attack  with  respect  to  the  non-rolling  coordinate  frame 

density 

coning  rate  of  projectile 

non-dimensional  coning  rate 

coning  rate  of  projectile,  non-dimensionalized  by 


Superscripts 

(•) 

(  )' 

(') 


rate  of  change  with  respect  to  time 

rate  of  change  with  respect  to  space 

quantity  is  referenced  to  the  non-rolling  coordinate  frame 


Subscripts 
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freestream  quantity 
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